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INCIDENT SHOCK-WAVE CHARACTERISTICS IN AIR, ARGON, 

CARBON DIOXIDE, AND HELIUM IN A SHOCK TUBE 
WITH UNHEATED HELIUM DRIVER 

Charles G. Miller IH and Jim J. Jones 
Langley Research Center 

SUMMARY 

Incident shock-wave velocities were measured in the Langley 6- inch expansion 
tube, operated as a shock tube, with air, argon, carbon dioxide, and helium as test gases. 
Unheated helium was used as the driver gas, and most data were obtained at pressures 

o 

of approximately 34 and 54 MN/m . A range of pressure ratio across the diaphragm 
was obtained by varying the quiescent test- gas pressure from 0.0276 to 34.5 kN/m^. 
Single- and double-diaphragm modes of operation were employed and diaphragms of 
various materials (aluminum, brass, copper, stainless steel, and plastic) tested. Shock 
velocities were determined from microwave interferometer measurements, response of 
pressure transducers positioned along the driven section (time -of -arrival gages), and 
measured tube-wall pressure. 

Comparison of the three methods for determining shock velocity showed that the 
microwave technique gives better definition of shock-wave characteristics (acceleration- 
deceleration trends) along the tube and is a viable method provided that microwave reflec- 
tion occurs at the shock front. However for the present range of flow conditions, the 
microwave technique proved to be rather limited in application. Shock velocities in- 
ferred from time -of -arrival transducers agreed within experimental accuracies with 
those inferred from measured wall pressures. The time-of-arrival method avoids many 
of the shortcomings of the other two methods and is believed to be the most acceptable 
method of those examined. Measured shock velocities were observed to increase ini- 
tially and then to attenuate with distance from the diaphragm for all test gases. This 
trend was most pronounced for helium; argon, air, and carbon dioxide experienced suc- 
cessively less attenuation. The distance required for the shock velocity to reach a maxi- 
mum value increased with increasing pressure ratio across the diaphragm for all test 
gases. At the higher pressure ratios, the shock velocity remained more nearly constant 
with distance after the initial acceleration for the present tube length of 29.0 m. Only 
minor differences in shock velocity were observed between single- and double-diaphragm 
modes of operation for a given test gas and pressure ratio across the diaphragm. The 



ratio of measured maximum shock velocity to predicted velocity was essentially the 
same for all diaphragm materials tested; however, the distance for the shock to reach 
a maximum value varied with diaphragm material. 

Measured maximum velocities exceeded predictions from conventional shock -tube 
theory for all test gases. Primary contributing factors to this discrepancy are believed 
to be (1) the effect of the helium driver temperature increasing upon pressurization of 
the driver section, (2) effects resulting from the finite opening time of the primary dia- 
phragm, and (3) two-dimensional or nonplanar flow effects. Tests with helium, which 
also yielded higher than theoretical shock speeds, suggest that interface mixing and flow 
chemistry are not primary contributors to this discrepancy. 

INTRODUCTION 

The Langley 6-inch expansion tube was designed for several modes of operation to 
provide a wide range of velocity-density performance. The first mode of driver opera- 
tion used in the facility after its initial operation made use of a resistance heater element 
mounted inside the driver chamber to heat the driver gas (refs. 1 and 2) to temperatures 
up to 600 K. Tests in this mode showed that incident shock speeds attained were in sat- 
isfactory agreement with theoretical one -dimensional prediction. At the conclusion of 
those tests, a series of tests was made without the driver-gas heater, and surprisingly 
the shock speeds achieved were only modestly lower than those with the heater. Also, it 
was noted that the tests without the heater showed significantly higher shock Mach num- 
bers than predicted by theory for a driver temperature of 300 K. 

As demonstrated in this report, the theoretical performance of the expansion tube, 
as predicted by conventional expansion- tube theory, is very sensitive to incident shock 
velocity in the intermediate section. Thus, in comparing predicted expansion-tube per- 
formance with measured performance, it is imperative that the input incident shock 
velocity be determined accurately. (It is well recognized that accurate measurement of 
the incident shock velocity in shock tubes is required for determination of the thermo- 
chemical conditions in the post- shock test gas.) 

The present investigation was undertaken to review and compare the accuracy of 
various means of determining incident shock velocity and to study more fully the phenom- 
enon of shock strengths which significantly exceeded theoretical prediction. For these 
tests the secondary diaphragm of the expansion tube was removed and the facility oper- 
ated as a conventional shock tube with an unheated driver gas. Air, argon, carbon dioxide, 
and helium were used as test gases because of their probable use in future expansion-tube 
studies. 


2 



Tests were conducted in two phases. After the initial testing phase, improvements 
were made to increase the precision of shock arrival time measurement, additional sta- 
tions were added, and a second testing phase was undertaken. The second phase covered 
a wider range of pressure ratio, included direct measurements of driver temperature, 
and since the facility configuration had changed in the interim, used a driven-section 
length of 29.0 m compared with 21.6 m for the first phase. 

Data from both series of tests are included in this report. Effects on shock- wave 
characteristics of (1) pressure ratio across the diaphragm, (2) test gas, (3) pressuriza- 
tion rate of the driver section (i.e., driver-gas temperature), (4) single- or double- 
diaphragm mode of operation, and (5) diaphragm material are examined. The various 
methods used to determine shock velocity are compared and relative advantages dis- 
cussed. Various nonideal effects are examined for their contribution to shock strength. 

SYMBOLS 

A ratio of cross-sectional area of driver section to cross-sectional area of 

driven section 

a speed of sound, m/sec 

Cp specific heat at constant pressure, J/kg-K 

d distance between successive minima or successive maxima on microwave 

film, m 

I distance from primary diaphragm station to position of maximum incident 

shock velocity, m 

Mg incident shock Mach number into quiescent test gas, Ug/aj^ 

Mg ^ incident shock Mach number calculated by compression model 

Ng electron number density, electrons/cm^ 

p pressure, N/m^ 

T temperature, K 
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t diaphragm thickness, m 

tg diaphragm thickness at groove, m (see fig. 3) 

U velocity, m/sec 

Ug incident shock velocity into quiescent test gas, m/sec 

V volume, m^ 

X axial distance downstream from the primary diaphragm, m 

n percent increase in shock Mach number (see eq. (4)) 

p density, kg/m 3 

T time, sec 

Subscripts: 

bf bottle field 

calc calculated 

d driver section 

i interface 

m measured 

max maximum 

op open 

1 state of quiescent test gas 

2 state of test gas behind incident shock in driven or intermediate section 

3 driver -gas conditions following unsteady expansion 
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driver- gas conditions at time of primary diaphragm rupture 


5 state of test- gas flow in acceleration section (free stream) 

5,t pitot pressure at expansion- tube test section 

Abbreviation: 

TOA time of arrival 

A bar over a symbol denotes average value. 

FACILITY AND APPARATUS 

The Langley 6-inch expansion tube is a cylindrical tube divided by two diaphragms 
(primary and secondary) into three sections. (See fig. 1.) The most upstream section 
is the driver or high-pressure section. This section is pressurized with a gas having a 
high speed of sound, such as helium. The intermediate section is evacuated and filled 
with the desired test gas at ambient temperature. The driver and intermediate sections 
are usually separated by a double-diaphragm apparatus, but a single diaphragm may also 
be used. The most downstream section is referred to as the expansion or acceleration 
section and has the same inside diameter (15.24 cm) as the intermediate section. A weak 
low-pressure diaphragm (secondary diaphragm) separates the driven and acceleration 
sections. The operating sequence of the expansion tube is shown schematically in fig- 
ure 1 and is discussed in reference 3. 

Results presented herein were obtained with the secondary diaphragm removed; 
that is, the expansion tube was operated as a shock tube. The length of the driver section 
was 2.44 m and the inside diameter was 16.51 cm. For the first series of tests, the 
length of the driven section was 21.6 m, whereas for the second series of tests, it was 
29.0 m. The inside diameter of the driven section was 15.24 cm, yielding a ratio of 
driver-section cross-sectional area to driven- section cross-sectional area of 1.174. 

Double-diaphragm mode of operation is generally employed to reduce the random- 
ness in pressure ratio across the primary diaphragm at the time of rupture. For this 
mode the driver section and double-diaphragm section (the volume between the two pri- 
mary diaphragms; see fig. 2) are pressurized with the driver gas to a pressure some- 
what less than the rupture pressure for a single diaphragm. The double -diaphragm 
section is then isolated from the driver section and the high-pressure supply field. The 
driver section is pressurized to the desired pressure, whereupon the double -diaphragm 
section is opened to atmospheric pressure. This decrease in pressure in the double- 
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diaphragm section results in rupture of the upstream diaphragm. Upon rupture of this 
diaphragm, the downstream diaphragm is subjected to a pressure on the driver-section 
side that is essentially that of the driver section. (The volume of the double -diaphragm 
section is small compared with the driver section, with the ratio of double- diaphragm - 
section volume to driver-section volume being approximately 0.07.) This pressure 
results in the rupture of the downstream diaphragm, and an incident shock wave is prop- 
agated into the quiescent test gas. 

Figure 3 shows the geometry for metal diaphragms. For most tests, diaphragms 
were fabricated from commercial plates of AISI type 304 stainless steel. A circular saw 
blade was used to cut cross-pattern grooves. In the present study, diaphragms used in 
the double -diaphragm mode were stainless steel and had a thickness t of 3.38 mm or 
4.775 mm, and the diaphragm thickness from the bottom of the groove to the side facing 
the driver section tg was 2.54 mm or 3.81 mm, respectively. For most tests made 
with a single diaphragm, t was 4.775 mm and tg was 3.81 mm. (For all single- 
diaphragm tests the diaphragm was placed in the downstream position in the double- 
diaphragm section (fig. 2), to maintain the same plane of initiation of shock formation 
(x equal to zero) as for the double -diaphragm tests.) Diaphragms with variations in t 
or tg exceeding 75 /J.m were rejected. 

Other metal diaphragm materials tested in the single -diaphragm mode of operation 
were aluminum, brass, and copper. The values of t, tg, and experimentally determined 
rupture pressure for the metal diaphragms are presented in the following table. Also 
presented in this table are the thickness and rupture pressure for a laminated plastic 
(polyethylene terephthalate) diaphragm. This diaphragm, which was unscribed, consisted 
of four layers of 0.36-mm-thick plastic and was positioned at the same station as the 
metal diaphragms. 


Diaphragm material 

t, 

mm 

tg, 

mm 

Rupture pressure, 
MN/m^ 

Stainless steel, AISI type 304 

4.78 

3.81 

34.9 

Stainless steel, AISI type 304 

3.38 

2.54 

20.7 

Stainless steel, AISI type 304 

2.08 

1.35 

9.7 

Brass, annealed 

3.53 

2.39 

6.7 

Brass, extra spring 

2.13 

1.27 

2.9 

Copper, cold roll, light 

4.45 

3.05 

5.4 

Aluminum, alloy 6061 

3.84 

2.57 

3.5 

Aluminum, alloy 2024 

4.93 

3.48 

5.4 

Polyethylene terephthalate, laminated 

1.44 


2.5 
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The square double -diaphragm section and the square-to-circular cross-section 
piece immediately downstream of the double -diaphragm section were modified to accept 
brass damper pads. (See fig. 2.) These pads have a thickness of 3.18 mm. The pads 
prevented the tips of the four diaphragm petals formed upon rupture from tearing off and 
thereby resulted in a significant reduction in debris carried downstream in the post-test 
flow. 


INSTRUMENTATION 
Incident Shock Velocity 


Microwave Interferometer System 

A microwave interferometer system (discussed in detail in ref. 4) was used to intro- 
duce microwave energy into the expansion tube, to excite a standing wave in the tube, and 
to record the movement of this standing wave due to the movement of the termination 
(region of ionized gas exceeding critical electron density) as a function of time. The 
varying periodic signal is recorded on a high-speed (60 rps) drum -camera-oscilloscope 
combination. Representative microwave films for argon, air, and carbon dioxide (CO2) 
are shown in figure 4. With helium as the test gas, the temperature immediately behind 
the incident shock was insufficient for ionization of helium. Hence, no microwave results 
were obtained, since the electron concentration was not high enough for reflection of the 
microwave signal. 

The antenna for this system for tests with a 21.6-m-long driven section consisted 
of an uninsulated tin-copper wire having a diameter of 0.61 mm. This wire was stretched 
horizontally across the tube diameter and positioned approximately 5 cm downstream of 
the tube exit. A carbon steel drill rod having a diameter of 4.77 mm and protruding 
5.08 cm into the tube was used as an antenna for tests with the 29.0-m-long driven section. 
This antenna was positioned 0.63 m upstream of the tube exit. The signal generator was 
adjusted to put out a 1.326-GHz signal. Film speed was checked for each test by utilizing 
a time mark generator to place marks on the film representing 100-|isec increments. 

Time-of-Arrival (TOA) Measurement 

A more conventional means of determining incident shock velocity Ug is to posi- 
tion high-frequency-response transducers along the length of the tube at known intervals. 
This procedure allows a distance-time history to be generated; hence, the average Ug 
is determined between successive instrumented stations. The time interval for incident 
shock arrival between stations was determined from the response of piezoelectric (quartz) 
pressure transducers mounted along the tube flush with the tube wall. At four stations. 
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thin-film resistance heat-transfer gages and photomultipliers with field stops to collimate 
the light source served as backup instrumentation. Station locations of pressure trans- 
ducers in terms of axial distance downstream from the most downstream diaphragm are 
presented in table I for both driven-section lengths. Outputs from these three types of 
instrumentation were recorded from an oscilloscope with the aid of a camera, and times 
for the shock to travel between stations were obtained from manual reading of the oscil- 
loscope film. Each type of instrumentation yields information about the flow behind the 
shock wave. With the exception of the pressure transducers, this instrumentation was 
used only to determine time of arrival of the shock wave and is referred to collectively 
as time -of -arrival transducers. For tests with the 29.0-m-long driven section, 

20 counter-timers were also used to obtain a time history of the incident shock, (See 
table I for locations of stations having counter-timers.) 

Reliable triggering of facility recording instrumentation is of primary concern. 

In the initial tests of the Langley 6-inch expansion tube, it was found that a reliable trig- 
gering signal for the oscilloscopes was produced by an accelerometer mounted to the 
exterior of the double-diaphragm section. This accelerometer detected an acceleration 
in the axial direction of the tube resulting from the diaphragm rupturing process. (For 
tests withihe plastic diaphragms the accelerometer did not provide an adequate trigger 
signal; thus the transducer output at station 8 (see table I) was used as a trigger signal.) 
For tests performed with the 21.6-m-long driven section, the signal from this acceler- 
ometer triggered oscilloscopes used to monitor the microwave signal and the output of 
pressure, heat-transfer, and photomultiplier transducers. Time delays to oscilloscopes 
monitoring transducers further downstream of the primary diaphragm were supplied by 
a digital time delay generator and are known to within 1 fxsec. This triggering method 
is not subject to spurious signals, nor does it require gain adjustments for varying pres- 
sure levels. 

To enable faster time sweeps of the oscilloscopes monitoring the transducer outputs 
along the entire tube for the second series of tests (29.0-m-long driven section), output 
signals from certain transducers were used as trigger signals. These signals triggered 
downstream oscilloscopes directly or through a digital time delay generator. This more 
complex triggering system, which required gain adjustments on charge amplifiers pro- 
cessing trigger signals from the pressure transducers for change in test conditions (vary- 
ing pressure level), performed with relatively good reliability for the present tests. The 
counter-timers used with the 29.0-m-long driven section were all started by the signal 
from the accelerometer mounted to the exterior of the double -diaphragm section, and 
each counter-timer was stopped upon shock arrival by the signal output from the pressure 
transducer at a given axial station. 
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Wall- Pressure Measurement 


The incident shock velocity at a given station may be computed from the measured 
wall pressure at that station. (The wall pressure is assumed to be an accurate represen- 
tation of the static pressure immediately behind the incident shock wave.) Thus, pressure 
transducers were also employed to measure the driven- section wall pressure with time. 
Each pressure transducer was used in conjunction with a charge amplifier and was cali- 
brated statically before the test series. The output signal was recorded from an oscil- 
loscope with a camera. Representative oscilloscope films illustrating wall-pressure 
variation with time at an axial station x of 19.75 m are presented in figure 5 for the 
three test gases employed in the first series of tests. Because of the larger number of 
pressure transducers used in the second series of tests and the fact that obtaining accu- 
rate pressure measurements is difficult and laborious, pressure instrumentation in the 
second series of tests was not treated with the same care. Hence, no attempt is made 
to infer incident shock velocity from wall-pressure measurements in tests with the 
29.0-m-long driven section. 

Driver -Section Pressure and Temperature 

Driver- section and double-diaphragm-section pressures were measured with strain- 
gage transducers statically calibrated to 68.95 MN/m2. Output from these transducers 
was read visually from a digital voltmeter and recorded on a strip chart. Quiescent 
test-gas pressdre was measured with a variable capacitance diaphragm -type transducer. 
Driver-section temperature was measured with a bare-wire chromel-alumel thermo- 
couple inserted through the upstream end plate of the driver section approximately 6 cm 
into the section. The thermocouple junction was exposed directly to the driver gas to 
provide the fast response required to obtain temperature histories during pressurization 
of the driver section. This thermocouple output was read from a compensated digital 
readout and recorded on a strip chart. Quiescent test-gas temperature was measured 
with a chromel-alumel thermocouple encased in a stainless-steel shroud. 

TESTS 

For the present tests the helium driver was drawn from a high-pressure supply 
field at ambient temperature, and no external heat was applied. The actual driver-gas 
temperatures obtained are discussed in a subsequent section. Approximate driver pres- 
sure at the time of diaphragm rupture p^ was 34 MN/m2 or 54 MN/m2 for the double- 
diaphragm tests and 35 MN/m2 for the single-diaphragm tests with stainless-steel 
diaphragms. Dry (bottled) air, and high purity argon, CO 2 , and helium were employed 
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as test gases. The range of pressure ratio across the primary diaphragm, for a given 
p^, was obtained by varying the qviiescent test-gas pressure Pj from 0.0276 to 
34.5 kN/m2. Tests performed with aluminum, brass, copper, and plastic diaphragms 
used helium as the test gas, and the pressure ratio across the diaphragm was approxi- 
mately 10^. The driven section was evacuated to approximately 0.01 N/m2 before filling 
with the test gas. 


DATA REDUCTION AND UNCERTAINTY 

Microwave Interferometer Measurements 

For the first series of tests with the 21.6-m-long driven section, the distance 
between successive minima or successive maxima d of the microwave traces was dig- 
itized with a film reader having a calibration factor of 30.11 micrometers per count. The 
repeatability in reading was within approximately ±4 counts, or ±120.44 micrometers. 

The range of velocity experienced yields d from roughly 3.05 mm to 5.08 mm; the cor- 
responding maximum uncertainty due to this error in reading between successive cycles 
is 4.0 and 2.4 percent, respectively. Although microwave traces obtained in the second 
series of tests were digitized with a film reader having a different calibration factor, 
the uncertainties due to film reading are expected to be the same as for the first series. 
However, because the microwave traces were read every other cycle in the second series 
of tests, the uncertainty in velocity due to reading of the film was reduced. By running 
a brass plug down an accurately determined length of the expansion tube and counting 
cycles, d was fovind to correspond to a movement of 22.71 cm. Since 5.08 cm of film 
corresponds to 1 msec, the time between successive minima or successive maxima is 
given by 


At = 


d 

5.08 cm/msec 


( 1 ) 


where d is expressed in centimeters. The corresponding velocity, in m/sec, is given 
by 


y _ 227.1 mm ^2) 

At 

where At is expressed in milliseconds. Uncertainties in the microwave system due 
to film speed variation, film resolution, and film length stability are discussed in refer- 
ence 4 and contribute a maximum error of less than 0.75 percent. Consideration of all 
sources of errors in measuring incremental velocities is believed to result in a maxi- 
mum uncertainty in velocity of 3 percent at the lower velocities of this study and 4.5 per- 
cent at the higher velocities. It should be emphasized that these uncertainties refer to 
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measurement of the velocity of an ionization front which has reached a critical electron 
density. (This critical density is 2.18 x 10^0 electrons/cm^ for the present microwave 
system .j An additional, possibly large, uncertainty may result if this critical density 
does not coincide with the incident shock front. 

Time of Arrival 

An average, velocity between stations was determined by knowing the distance 
between stations and the time for the shock wave to travel this distance. These times 
were obtained from the response of pressure transducers along the driven section, the 
output of which was recorded on oscilloscope films and displayed on counter-timers. 
Uncertainties in pressure transducer response, oscilloscope time scale, and reading of 
oscilloscope film are believed to result in a corresponding uncertainty in time interval 
between arrival of the shock wave at successive stations of less than ±10 psec for the 
first series of tests and ±5 psec for the second series. Uncertainty in oscilloscope time 
scale was reduced by using a timing mark generator to supply a known time increment 
to the oscilloscope. If necessary, a correction was applied to film readings. With one 
exception, oscilloscope time scales were accurate to within 3.75 percent. The principal 
sources of error for times obtained from the counter-timers are response of the pres- 
sure transducers and associated equipment and counter-timer sensitivity to this response. 
For example, at the lower values of incident shock velocity, disturbances resulting from 
rupture of the thick steel diaphragms propagated through the tube wall and preceded the 
arrival of the incident shock. In this case, many of the pressure transducers responded 
to these acceleration effects, thereby producing a premature stop of the counter-timers. 

At the higher values of incident shock velocity, the static pressure level behind the inci- 
dent shock was often too low to provide sufficient output from the charge amplifier of the 
pressure transducer to stop the counter-timer. In this case the stop would correspond 
to flow in the higher pressure region of the expansion fan and hence yield a delayed stop 
signal. Counter times which were obviously erroneous were eliminated from the data, 
and the velocity was averaged between the adjacent stations. For both the oscilloscope 
films and the counter-timer readings, uncertainties in shock velocity inferred from TOA 
transducers depend on the distance between successive stations, being larger for stations 
closer together. The zero time for these distance-time plots was also obtained from the 
accelerometer signal. Randomness involved in the accelerometer signal results in a 
shift in the time scale between individual tests. 

Wall Pressure 

A third method employed to infer Ug was to use the measured driven-section 
wall pressure P 2 as input to the shock-tube phase of the computer program of refer- 
ence 5 (in conjimction with measured p]^ and Tj). Experimental uncertainties in P 2 
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are dependent on many factors, such as transducer sensitivity, calibration technique 
(static or dynamic), transducer linearity, transducer mounting, adequate thermal protec- 
tion, oscilloscope accvtracy, quality of oscilloscope traces with respect to the signal- 
noise ratio, and oscilloscope film reading procedure. The maximum uncertainty in P 2 
for the tests in the first series is believed to be less than 10 percent. For the conditions 
of this study the uncertainty in Ug is roughly half that of the corresponding P 2 - 

RESULTS AND DISCUSSION 

Effect of Uncertainty in Incident Shock Velocity on Predicted 
Expansion- Tube Performance 

For any shock tube or shock -tube -related device such as the expansion tube, it is 
well-known that an accurate determination of the incident shock velocity Ug is required 
for determination of the thermochemical state of the post-shock gas. In expansion -tube 
operation, measurement accuracy of Ug acquires additional importance because of the 
multiplying effect of the subsequent expansion. The expansion fan which passes through 
the shocked gas raises the gas to a higher velocity while expanding it to a lower pres- 
sure and temperature. This final state after expansion is dependent on the enthalpy of 
the gas entering the fan. Rather small differences in this initial enthalpy can result in 
large differences in the final state of the gas exiting the fan. 

Examples of the effect of uncertainty in Ug on predicted expansion -tube flow 
quantities for air, argon, CO 2 , and helium as test gases are shown in figure 6. The free- 
stream static pressure P 5 and pitot pressure P 5 ^ are shown as functions of free- 
stream velocity U 5 . These results were calculated with the conventional expansion -tube 
program of reference 5 for thermochemical -equilibrium flow and no shock reflection at 
the secondary diaphragm. The calculations are for values of incident shock velocity Ug 
which correspond to those attained in the Langley 6 -inch expansion tube for an unheated 
helium driver pressure of 34 MN/m^. An uncertainty in Ug of 5 percent is assigned 
for each gas. 

The effect of this uncertainty in Ug on predicted free -stream static pressure and 
pitot pressure is illustrated by the shaded region in figure 6 . The region of uncertainty 
is magnified as the gas expands to higher velocity. At a velocity of 5.5 km/sec, typical 
of expansion -tube operation, the uncertainty in static pressure ranges from about 30 to 
70 percent for the various gases, and the uncertainty in pitot pressure varies from about 
20 to 45 percent. Thus, relatively small uncertainties in input Ug result in much 
larger uncertainties in predicted P 5 and pg ^ for a given flow velocity. Hence, to 
minimize uncertainties in predicted expansion -tube test-section flow conditions, the inci- 
dent shock velocity Ug should be measured as accurately as possible. 
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II 


Comparison of Methods for Measur ing Incident Shock Velocity 

TOA data from counter -timer readings and from oscilloscope film readings are 
presented in figures 7, 8, 9, and 10 for air, argon, CO 2 , and helium as test gases, respec- 
tively. In these figures the incident shock velocity Us is plotted as a function of axial 
distance x downstream from the primary diaphragm for various values of quiescent 
test-gas pressure pj. Horizontal filled bars denote the average velocity between suc- 
cessive stations from oscilloscope film; horizontal open bars denote the average velocity 
determined from counter -timer readings. Also shown in these figures are values of Ug 
determined from microwave measurements. These microwave data were read every 
other cycle and are connected by line segments. The results of figures 7 to 10 were 
obtained with the double -diaphragm mode of operation and a helium driver pressure of 
approximately 34 MN/m^. 

Microwave results were not obtained at the higher values of quiescent test-gas 
pressure pj for air (fig. 7) and CO 2 (fig. 9), nor for helium (fig. 10) at any test value 
of pj. This absence of microwave data at these conditions is due to Insufficient electron 
density at the shock front or in the post -shock flow region to terminate the microwave 
signal. Ideally, TOA data from counter -timer readings should agree with data from 
oscilloscope films for a given test. The differences observed in figures 7 to 10 between 
these two TOA methods are indicative of the present experimental uncertainties, which 
in turn are dependent on limitations of the methods. For example, at the higher values 
of pj for all test gases (corresponding to the lower values of Ug), disturbance created 
by the rupture of the thick steel diaphragms propagated downstream through the tube 
wall and preceded the arrival of the incident shock. Several pressure transducers 
experienced rather severe acceleration effects because of this propagation, and the 
resulting output from the transducers caused a premature stop of the counter -timers. 
These same acceleration effects on pressure transducer outputs used to trigger down- 
stream oscilloscopes resulted in erroneous triggering. The high noise -signal ratio also 
led to difficulty in reading the time of arrival of the incident shock from oscilloscope 
film. At the lower values of pj, the incident shock preceded disturbances in the tube 
wall, but the pressure immediately behind the shock was sometimes insufficient to pro- 
vide adequate output from the charge amplifier of the pressure transducer to stop the 
counter -timers on arrival of the shock or to trigger oscilloscopes. 

Relatively good quality microwave film traces were obtained for air at values of pj 
less than 3.45 kN/m^. As observed in figure 7, the microwave method provides more 
detailed information on the shock travel, since an amplitude peak occurs for every 
2 2. 71 -cm movement of the front being tracked. This small increment, which is equal 
to 45.42 cm in figures 7 to 9 since every other cycle of the microwave trace was read, 
means that a given uncertainty in the time increment At results in relatively large 
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uncertainties in the corresponding velocities. The percent uncertainty in At increases 
with increasing shock velocity (see the section entitled "Data Reduction and Uncertainty"); 
hence the scatter in microwave data is expected to increase with decreasing pj^. At 
values of pj for which good microwave traces were obtained in air, the shock velocities 
from the microwave measurements generally were somewhat lower than shock velocities 
obtained from the TOA method. 

Microwave film traces for the argon tests displayed good signal amplitude for all 
values of quiescent test -gas pressure. As shown in figure 8, shock velocities from the 
microwave method and TOA method are in good agreement for values of pj between 34.5 
and 3.45 kN/m^. However, at the higher shock velocities corresponding to values of pj^ 
of 0.69 and 0.345 kN/m2, large discrepancies exist between the two methods for deter- 
mining Ug. The microwave results for these two values of pj^ are suspect, from the 
viewpoint of the trend of Ug with x and the fact that the peak values of Ug exceed 
the theoretical limiting value (the value of Ug as P4/P1 approaches infinity^. As the 
quiescent pressure is decreased further to 0.069 and 0.0345 kN/m^, the agreement 
between the two methods improves markedly, but the scatter in Ug inferred from 
microwave measurements is significantly larger than the scatter observed at the higher 
values of pj, as expected. 

Microwave film traces for CO2 tests were characterized by low amplitude and were 
frequently irregular. As for air, shock velocities from the microwave method for CO2 
(fig. 9) were generally less than those determined from TOA data. The absence of micro- 
wave data for helium (fig. 10) is expected, since (1) the range of shock velocity was insuf- 
ficient to ionize the helium and (2) the photomultipliers failed to show any light emitted 
from the helium, and thus the test gas was assumed to be essentially free of contaminants 
which are more readily ionized. 

TOA data and microwave data can also be conveniently compared on a plot of time 
of incident shock arrival as a function of axial distance downstream of the primary dia- 
phragm. Such plots, referred to herein as distance -time diagrams, are shown in fig- 
ure 11 for air, argon, and CO2 and for selected values of quiescent test-gas pressure. 

The microwave data were forced to agree with the TOA data at the 5.377-m station or 
11.640-m station. (Because the accelerometer was used as a trigger for the microwave 
interferometer system, the microwave data lack an absolute reference point. On the 
two drum -camera-oscilloscope systems used to record the microwave signal, the signal 
from the pressure transducer at x equal to 5.377 m was placed on one film and the sig- 
nal from the pressure transducer at x equal to 11.640 m on the other.) Although the 
scales of figure 11 are relatively coarse (for example, the symbol width for TOA data 
corresponds to approximately 2 cycles on the microwave trace), these distance -time dia- 
grams show that TOA and microwave measurements do not generally agree as to the 
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position of the shock front at a given time over the entire distance traversed by the 
shock. Examples of particularly poor agreement are shown for argon at pj equal to 
0.345 kN/m2 (fig. 11(b)) and for CO 2 at this same value of pj (fig. 11(c)). The observed 
discrepancy between microwave and TOA data at other values of p^ for these two test 
gases and for air exceeds the uncertainties in measurement. Although the point of forced 
agreement on the distance -time diagrams is arbitrary, selection of a different adjustment 
point will not bring the data into agreement for the entire traverse of the shock wave. It 
seems obvious that the microwave signal sometimes reflects from a region other than the 
pressure front. 

As noted previously, microwave measurements with the present system were not 
possible for helium tests and tests with air and CO 2 at the higher values of quiescent 
test-gas pressure. In table II, the maximum and final (near driven-section exit) inci- 
dent shock velocities for a range of pj for each test gas are tabulated. These results 
were obtained from the TOA data of figures 7 to 10. Also tabulated are predicted elec- 
tron number densities corresponding to the maximum and final shock velocities for a 
given pj. These electron densities were calculated by the thermochemical -equilibrium 
program of reference 5. For the microwave signal to be terminated, an electron density 
of 2.18 X 10^0 electrons/cm^ is required. Hence, the results of this table illustrate that 
insufficient electron densities were present to provide microwave termination for the 
cases previously mentioned. According to the data of table II, sufficient electron density 
existed for the argon tests. However, at values of p^ of 0.69 and 0.345 kN/m^, the 
microwave data exhibited a departure from the TOA data not observed for the other test 
gases, nor for other values of pj^ for argon. At these two values of pj, the microwave 
signal is terminated at something other than the shock front. 

Immediately behind the shock front, translational degrees of freedom are excited 
within several collision times, whereas excitation of internal degrees of freedom, dissoci- 
ation, and ionization take many more collisions. Thus, a relaxation zone is produced 
immediately behind the shock front. For argon, electron production proceeds initially by 
atom-atom collisions and collisions of electrons from impurities with atoms until suffi- 
cient electrons exist for electron-atom collisions. Hence, relatively few electrons exist 
at the shock front, and the equilibrium degree of ionization can occur at a significant dis- 
tance behind the shock front. (See, for example, refs. 6, 7, and 8.) Measvmements 
associated with the chemistry of shocks into argon are beyond the scope of this study. A 
possible explanation for the difference observed between the microwave and TOA data 
for argon at pj of 0.69 and 0.345 kN/m^ is that these conditions correspond to an 
appreciable departure from equilibrium resulting in an electron density deficiency in the 
relaxation zone. At the lower values of pj^, the role of impurities in generating elec- 
trons may become more significant. (See refs. 7, 8, and 9.) The supply of argon used 
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in the present study was 99.998 percent pvire, but outgasing of the tube and dump tank 
walls and leaks in the vacuum system may have introduced a significant impurity level. 
This level could have an appreciable effect on the microwave measurements in argon, 
where high temperatures were generated in the post -shock flow region. (See table II. ) 
Whereas flow impurities may have a significant effect on the electron density, the corre- 
sponding effect on shock velocity and thermodynamic properties in region (see fig. 1) 
may be negligible. For example, a test relating to another investigation was performed 
in which a mixture of 98 percent argon and 2 percent methane was used at a pj of 
6.9 kN/m^ and p^ of 34 MN/m^. With the exception of the photomultiplier traces, no 
appreciable difference was detected between this test and a similar test with pure argon. 

Prior to a series of expansion-tube tests, pressure transducers mounted along 
the tube wall are usually dynamically calibrated by employing the facility as a low- 
performance shock tube. For such checks on the static calibration of the transducers 
and on thermal insulation of the transducer sensing surface, a plastic diaphragm is used 
to avoid generation of tube -wall disturbances and to minimize cost of operation. This 
procedure was reversed in the present study in that shock velocities were determined 
from measimed wall pressures. For a given axial station, the pressure measured 
immediately behind the incident shock was used as input to the program of reference 5 
to calculate the corresponding shock velocity. In figure 12, values of Ug calculated 
from tube -wall pressure measurements for air, argon, and CO 2 are compared with shock 
velocities determined from microwave measvirements and TOA data. These sample data 
were obtained with the 21.6-m-long driven section, double -diaphragm mode of operation, 
and a helium driver pressure of approximately 34 MN/m^. For the most part, values of 
shock velocity inferred from TOA data and pressure data agree within the experimental 
accuracy, as do the microwave data for air and CO 2 . 

Thus, comparison of the three methods for determining shock velocity shows that 
the dependence of the microwave technique on electron density at the shock front limits 
its range of applicability for a study such as this one employing a number of test gases 
over a range of shock velocity. Use of the microwave technique should always be accom- 
panied by TOA transducer information to verify its validity, since good quality microwave 
film traces (indicating termination of the microwave signal) did not always yield a posi- 
tion of reflection that was in agreement with the pressure front. In its range of validity, 
the microwave technique is a useful tool because it yields better definition of shock -wave 
acceleration-deceleration characteristics than practical with TOA transducers. Of the 
two other methods, the TOA method is believed to be more desirable than the pressure 
method, since less effort is required to obtain shock velocity with this method and the 
pressvire method has less potential for improvement. Factors contributing to the greater 
effort for the pressure method are: (1) frequent and careful calibrations of the pressure 
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transducers and associated charge amplifiers and oscilloscopes are required; (2) care 
must be exercised in moxmting the transducers to insure that the sensing surface is flush 
with the tube wall, that adequate thermal insulation for the transducers is provided, and 
that the mounting procedure reduces acceleration effects resulting from disturbances 
transmitted along the tube wall; and (3) signal noise levels and trace characteristics 
require judgments to be made in reading the outputs. The inferred pressures must then 
be used as inputs to a program such as reference 5 to determine corresponding values of 
shock velocity. The accuracy of this technique is believed to be limited to about 5 per- 
cent in shock velocity. 

Although pressure transducers were used herein for the TOA method, thin-film 
resistance transducers are preferable, since they have a faster response to shock arrival 
and a wide range of flow conditions for which they give good definition of the shock posi- 
tion, they are not influenced by disturbances propagating through the tube wall, and they 
provide a better signal -noise ratio. However, wall -pressure measurements are basic in 
expansion-tube testing, and it was not feasible to replace existing pressure transducers 
with heat -transfer transducers in the Langley 6 -inch expansion tube for the present 
shock -tube tests. 


Comparison of Measured Shock Velocity With Theory 

As observed for all four test gases (figs. 7 to 10) and in most shock -tube studies 
(for example, refs. 10 to 14), the incident shock initially accelerates for some distance 
down the tube, reaches a maximum value Ug uiax’ then decelerates for the remainder 
of its travel. The principal cause of the acceleration phase is the finite opening time of 
the diaphragm, which causes the shock wave to be formed and strengthened by a family of 
compression waves which emanate from the diaphragm location as the open area increases 
with time. The dominant cause of the deceleration phase is the growth of the wall bound- 
ary layer, although the diaphragm opening process may also make a small contribution to 
this deceleration. 

As the higher values of pressvure ratio across the diaphragm P4/P1 are approached, 
the incident shock velocity tends to remain more constant with distance x after the 
initial acceleration for the driven-section length of 29.0 m. The deceleration phase is not 
examined in detail in the present study. Linear fits to the downstream portion of the 
velocity -distance data of figures 7 to 10 revealed that helium suffers the most severe 
attenuation. Argon, air, and CO2 experience successively less attenuation. 

The distance I required for the shock wave to reach Ug is a fimction of the 

time required for the diaphragm to open. The characteristic wave which originates from 
the diaphragm location at the instant of complete opening may be viewed as the completion 
of the shock-strengthening wave family. The distance downstream of the diaphr^m 
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where this characteristic wave overtakes the incident shock wave represents the dis- 
tance I, provided attenuation effects are ignored. An approximate method of estimating 
the distance I (for example, refs. 11 and 12 ) has been employed in which l is assumed 
to be proportional to the product of the diaphragm opening time and the maximum shock 
velocity. The constant of proportionality is determined empirically and the resulting 
simple relations from references 11 and 12 are shown graphically in figure 13. The 
shaded region of figure 13 denotes the present data for all four test gases. The relatively 
large size of the shaded region is due to data scatter and imcertainty in defining I and 
Us from figures 7 to 10. Diaphragm opening time T^p was not measured in the 

present study, and predicted values for the stainless -steel diaphragms varied consider- 
ably according to the method of calculation. For the present diaphragm thickness and 
opening area, the method of reference 15 is believed to provide a reasonably accurate 
value of Tpp. The diaphragm opening time predicted from the method of reference 15 
is approximately 500 iisec, and this value is assumed in presenting the data of figure 13. 
The dashed lines of figure 13 denote the regions for which the empirical relations of ref- 
erences 11 and 12 were obtained. As observed from this figure, the present results 
represent a large extension in the range of I of these previous studies. At the lower 
values of pressure ratio P 4 /P 1 (that is, lower values of Ug) of this investigation, the 
present data and the relation of reference 11 are in reasonably good agreement. How- 
ever, as the pressure ratio increases, the present values of I increase and for a given 
^op^s max’ ®xceed values calculated from the relation of reference 11. Because of the 
uncertainty in predicted opening time TQp, as well as uncertainties in determining 
Ug jnax definite conclusions may be drawn. However, this comparison sug- 

gests that the relations of references 11 and 12 are not valid over the wide range of test 
conditions obtained in the present facility and should be used with discretion. 

Comparison With Conventional Theory 

The dashed lines in figiores 7 to 10 represent the predicted shock velocity from 
conventional shock-tube theory for each given diaphragm pressure ratio. This conven- 
tional theory assumes one -dimensional, constant -area, inviscid flow and thus includes no 
mechanism for shock attenuation. Instantaneous diaphragm rupture and shock formation 
are assumed, and both driver -gas and driven-gas initial temperatures (T^ and Tj, 
respectively) were taken to be 298 K. The computations assumed equilibrium real -gas 
properties for the shocked test gases and considered imperfect -gas effects for the helium 
driver. In all cases, the measured maximum shock velocity exceeds prediction. Because 
of attenuation, which is attributed primarily to viscous effects, agreement between pre- 
dicted and measured velocities improves near the tube exit, but in many cases the final 
velocity at the end of the tube still exceeded the prediction from conventional theory. 
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Experimental shock velocities which exceed conventional shock -tube theory have 
been previously noted in the literature. For example, in reference 16 for diaphragm 
pressure ratios larger than 10®, the measured shock velocities obtained with a combus- 
tion driver approached twice the theoretical value. In reference 17 with an unheated 
helium driver, measured velocities exceeded conventional shock -tube theory by as much 
as 25 percent. In reference 10 unheated helium and hydrogen driver gases are shown to 
produce higher shock velocities than predicted for nitrogen, oxygen, air, or argon as 
driven gases. 

A far greater number of sources could be cited for which the shock velocity does 
not exceed prediction, so apparently the conditions which cause this effect are not usually 
dominant. In the following section, several possible causes for the shock velocity being 
greater than predicted are examined. 

Causes of Discrepancy 

Driver temperature . - It is well-known that the driver temperature strongly 

influences the shock strength. (See ref. 18.) Since both the high-pressure helium stor- 
age bottle field and the driver section were at ambient temperature, the ambient tempera- 
ture of the driver section was used in the predictions illustrated in figimes 7 to 10. 
However, the Joule -Thomson coefficient is negative for helium in this temperattu'e range, 
and the expansion of helium from the storage field into the driver section results in a 
temperature rise. In addition to the Joule -Thomson heating, a compression heating 
occurs since the helium in the driver section is continuously being compressed during 
the charging process. If this additional temperature is not lost by heat transfer to the 
driver-section walls prior to diaphragm rupture, a more effective driver condition 
results. The temperature -time history of the helium during pressurization of the driver 
section was monitored for each test with the 29.0-m-long driven section. In figure 14 
the maximum driver temperature recorded during pressurization is plotted as a function 
of supply pressure in the helium bottle field for a number of tests. This figure demon- 
strates that values of T^^ during pressurization were frequently over 340 K. A 

least-squares fit to the data of figure 14 illustrates the trend of increasing T^j j^ax '^ith 
increasing bottle -field pressure. 

An example of driver -section pressure and temperature variation with pressuriza- 
tion time is shown in figure 15. For this example, the driver section was charged to 
about 30 MN/m^ and then isolated from the bottle field for an additional 2 min. During 
the charging time the driver -gas temperature, as measured by the thermocouple, 
remained between 330 K and 340 K. When the supply valve to the bottle field is closed, 
the temperature decays to the ambient driver -section wall temperature because of heat 
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loss. Since the density is constant after the valve is closed, a corresponding decrease 
in pressure occurs. 

In normal double -diaphragm operation the diaphragm rupture generally occurs 
approximately 5 to 10 sec after closing the supply valve. On each rim, driver pressure 
and temperature were noted at the time of diaphragm rupture. These values of pres- 
sure p^ and temperatvme T4 were used as inputs in the theoretical computations of 
shock velocity by the method of reference 5 and are indicated in figures 7 to 10 by the 
long-short dashed lines = T4 m)* 

The calculated effect of driver temperature for all four test gases is shown in fig- 
ure 16, in which the ratio of driver -section pressure to quiescent test-gas pressure is 
plotted as a function of incident shock Mach number Mg = Ug^a^. The measured values 
of T4 and P4 were used as input to the conventional shock -tube method of reference 5 
to predict the shock Mach number for each test. The solid lines of figiore 16 represent a 
fairing of these predicted values of Mg. For comparison, similar predictions using a 
value of T4 of 298 K are shown, along with the measured maximum shock Mach num- 
ber. The effect of driver temperature is observed to be insxifficient to explain the dif- 
ferences between measured and predicted shock velocities. Also shown in figure 16 are 
the measured maximum shock Mach numbers for double -diaphragm tests in which the 
driver pressure was increased to approximately 54 MN/m^. These higher values of P4 
were obtained to extend the range of P4/P1 and to determine whether any effect on max- 
imum shock Mach number resulted from operation at a higher driver pressure. No sig- 
nificant effect of pressure level on shock velocity is observed. 

It is interesting to note that measured values of maximum shock Mach number for 
air (fig. 16(a)) and CO2 (fig. 16(c)) are amenable to a straight-line fit on a semilogarith- 
mic plot. The resulting expressions are: 

For air. 


“s,max “ 0-92 ^ 


(3a) 


For CO2, 


M 


s,max 


Pa 

« 1.2 logg ~ - 1.56 


(3b) 


These relations for air and CO2 describe the experimental results presented in figure 16 
to within 3.5 and 3 percent, respectively. Such simple expressions may be readily 
incorporated into the program of reference 5, and thereby permit a more accurate pre- 
diction of expansion -tube performance for these two test gases than previously possible 
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by using inputs and to determine the incident shock velocity in the intermediate 
section. It should be emphasized that these relations are not general and discretion must 
be exercised in usage for other shock tubes. 

The effect of driver temperature on shock velocity was verified experimentally by 
varying the pressurization time of the driver section. A number of tests were performed 
in which the driver gas was allowed to return to near ambient temperatvire before initiat- 
ing diaphragm rupture. For these "slow-charge" tests, additional helium was added as 
needed to make up the pressure loss due to cooling. The resultant shock velocities from 
these slow-charge tests are compared in figvure 17 with "fast -charge" tests with air and 
helium and for several values of quiescent test -gas pressure. These data illustrate that 
the higher driver temperature due to fast charge increases the shock velocity, as expected. 
The velocity difference between the two charge rates is essentially constant down the 
length of the driven section, and shock attenuation does not appear to be significantly 
affected by this difference in driver temperature. As observed from figure 7 for air and 
figure 10 for helium, conventional shock -tube theory predicts a difference in shock veloc- 
ity (Mach number) roughly the same as observed experimentally for this temperature 
change. Although elevated driver temperature is a contributing factor, it does not fully 
account for the discrepancy between experiment and theory. 

Single- or double -diaphr agm mode of operation .- With the double -diaphragm mode 
of operation, the test is initiated by venting the chamber between diaphragms, which 
causes first the upstream and then the downstream diaphragm to rupture, since each one 
successively feels the entire driver pressure. This arrangement is not unlike a three - 
chamber shock tube which is generally called a bviffered shock tube. Buffered shock 
tubes have been shown (ref. 19) to produce stronger incident shocks in the driven section 
than single -diaphragm shock tubes for equivalent driver- and driven-section initial con- 
ditions. This effect seems unlikely in the present tests since the chamber between the 
diaphragms is only 15 cm long. In order to check this effect, a number of tests were 
made for all fovir test gases with a single diaphragm. These results are compared in 
figure 18 with double -diaphragm results. For several tests shown in figure 18, particu- 
larly at the highest diaphragm pressure ratios, the shock velocities for double -diaphragm 
tests reach a higher peak somewhat sooner (that is, at slightly smaller distance traveled). 
However, as a group, the shock velocities are in good agreement for the two operation 
modes. It is concluded that the diaphragm mode used has only minor effects on shock 
velocity. 

To learn whether the driver expansion process was significantly affected by the 
double -diaphragm arrangement, a pressure transducer was mounted in the driver -section 
wall a short distance (12.7 cm) upstream of the most upstream diaphragm station. Fig- 
ure 19 compares the pressure histories for single -diaphragm and double -diaphragm 
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operation. The data were obtained from a digital data storage unit operated in the 
’’pretrigger" mode, wherein the stored data include those measured just prior to the 
trigger signal. This allowed the preservation of the base line and early pressure decay 
even thoi^h the trigger signal was derived from the pressure decay. 

The dotted curve in figure 19 is the experimental driver pressure for a single - 
diaphragm test. Time equal to zero is an arbitrary time prior to diaphragm rupture. 
The pressure decay down to about 15 MN/m^ is observed to be in good agreement with 
that calculated for a centered unsteady expansion (dashed curve). At about 15 MN/m^ 
the pressure levels out, since, because of the larger cross-sectional area of the driver 
compared with the driven section, the unsteady expansion is not continuous to a Mach 
number of 1. (See the following sketch.) The pressure begins to decay again when the 
expansion, reflected from the upstream end of the driver section, crosses the pressure 
port. 



The solid curve of figure 19, taken from a double -diaphragm test, is in good agree- 
ment with the single -diaphragm results if the time scales are alined as shown. The first 
pressure drop with the double diaphragm undoubtedly occurs as a result of the rupture of 
the upstream diaphragm. A recompression occurs because the second diaphragm has 
not yet ruptured; but when it ruptures, the expansion process is very similar to that for 
the single diaphragm. The weak shock which moves upstream as a result of the recom- 
pression has apparently been attenuated before passing the pressure port again. 

These results illustrate the similarity of the driver expansion process for single - 
and double -diaphragm operation and support the previous conclusion that the diaphragm 
mode of operation has little effect on shock strength. 

Diaphragm material . - Several tests were performed with single diaphragms of dif- 
ferent material and thickness. Results of these tests are shown in figm-e 20, in which the 
average incident shock velocity between successive stations is plotted as a function of 
distance downstream of the diaphragm. Values of t, tg, and experimentally determined 
rupture pressure for the various diaphragms are presented in a previous section entitled 
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"Facility and Apparatus.” Whereas the scribed aluminum, brass, copper, and stainless - 
steel diaphragms opened into and conformed with the square -to -round transition section 
shown in figure 2, the unscribed plastic diaphragms opened into an essentially round ori- 
fice and thus did not conform to the wall of the transition section. Helium was employed 
as the test gas because of its ideal -gas behavior and because of the fact that of the four 
test gases examined, only helium exhibited significant differences between measured and 
predicted values of maximum incident shock Mach number at pressure ratios as low 

O 

as 10 . (See fig. 16.) Because of the lower rupture pressure of the diaphragms con- 
sidered in figure 20, it was necessary to test at a low pressure ratio (hence relatively 
high value of pj) to avoid instrumentation difficulties with the TOA method as discussed 
previously. For all tests presented in figure 20, the driver temperature was within 
2 percent of 333 K and the pressure ratio P4/P1 was within 9 percent of 1030. The 
initial acceleration to a maximum shock velocity and subsequent deceleration are observed 
for all diaphragm materials. The stainless -steel diaphragm and the plastic diaphragm 
yield the maximum shock velocity at a smaller distance I downstream of the diaphragm 
than the other three diaphragm materials. The expression of reference 11 was used to 
predict estimates of relative diaphragm opening times for the various materials. The 
fact that diaphragms having the fastest predicted opening times (plastic and stainless steel) 
yield the smallest values of I indicates an effect of diaphragm opening characteristics. 
Comparison of experiment with prediction (dashed lines of fig. 20 represent conventional 
shock -tube theory with measured p^, T4, pj, and Tj^ Inputs) reveals that the mea- 

sured maximum shock velocity is 1.08 to 1.11 times the predicted shock velocity for all 
diaphragms. The lowest values of Ug max/^s calc with brass and copper dia- 

phragms. For all diaphragms, the measured shock velocity exceeds prediction in the 
first 16 to 22 m of travel. 

Diaphragm opening time . - As noted previously. White (ref. 10) reported shock 
velocities which were greater than predicted by conventional theory. White attributed 
this difference to the influence of the finite time required to open the diaphragm. Con- 
ventional theory assumed that the diaphragm opens instantly, with Immediate shock for- 
mation. However, if a diaphragm is relatively thick and massive, it may require several 
hundred microseconds to open fully. During this opening, compression waves are gen- 
erated which coalesce into a shock some distance downstream. The resulting shock wave 
is stronger than the instantly formed shock. 

Formation -from -compression predictions: White solved the special problem in 
which all the compression waves formed a centered compression fan. The interface was 
assumed to accelerate in such a way that the compressions all coalesced at the same 
axial station to form a shock wave. This assumption simplifies the calculation since the 
driven gas is processed either by Isentropic waves only or by shock only. 
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Ikui, Matsuo, and Nagai (ref. 20) solved a more general problem, in which an arbi- 
trary number of convergence points were allowed. The accompanying sketch illustrates, 
in the distance -time plane, these two solution methods. 

I Interface 
s Shock wave 
c Compression wave 
r Reflected wave 

Centered compression (ref. lO) Multiple-step compression (ref. 20) 




The following assumptions inhibit the accuracy or completeness of these solution 
methods: 

(1) The reflection waves, denoted r in the sketch, which result from the inter- 
action of the compression wave and the shock, are computed, but the interaction of this 
reflection with the interface is ignored. This interaction results in a change in the inter- 
face speed, a transmitted upstream -facing wave into the driver gas, and a downstream - 
facing wave which catches the shock and modifies its strength. 

(2) The interface is assumed to be the generating surface of the compression 
waves, whereas in fact the waves originate at the diaphragm station as a result of varia- 
tion in area as the diaphragm opens. This assumption causes negligible error in the cal- 
culation of shock strength, but it does prevent a description of the wave process in the 
expanded driver gas and divorces the interface acceleration from its cause, the diaphragm 
opening. 

(3) White (ref. 10) uses real -gas properties for the driven gas, but Ikui, Matsuo, 
and Nagai (ref. 20) confined their calcxalations to ideal gases. None of these authors per- 
formed their calculations for CO 2 or helium as driven gases. 
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The function II was defined in reference 20 as the percent increase in shock 
Mach number due to formation from compression. Figure 21 shows the values of II for 
driver -driven gas combinations of helium -air, helium -argon, helium -CO2, and helium - 
helium. Both the method of reference 10 (single compression step) and that of refer- 
ence 20 (multiple compression step) are shown in figure 21. The curves for helium-air 
and helium -argon are taken from reference 20 and those for helium -CO2 and helium - 
helium were computed. Since II is defined by 


n s 




Me 



(4) 


the corrected shock Mach number at a given diaphragm pressure ratio P4/P1 is given 
by 


^s,c ~ ^s 



Although derived for A = 1, T4 = 300 K, and perfect gases, the II correction factors 
have been applied in the present case for A = 1.17, T4 = 330 K, and real gases. (The 
effects of an area change near the diaphragm are discussed in ref. 21.) These various 
predictions of shock Mach number as a function of diaphragm pressure ratio are shown 
in figixre 22, together with the measured values of maximum shock Mach number from 
both phases of the present study, double- and single -diaphragm modes of operation, fast 
and slow charge of driver section, and two values of driver pressure. Although mea- 
sured values of Mg are in better agreement with predictions from the formation- 

from -compression theories, the experimental values of Mg ^^ax’ the exception of 
the argon data, still exceed these predictions at the higher values of P4/Pl* 

Unsteady flow analyses: The method of characteristics was used to obtain a solu- 

tion for the flow development with a finite diaphragm opening time. The purpose was to 
compare the predicted location of the maximum shock velocity with experiment and to 
compare the measured pressure rise in the driver expansion fan with theoretical predic- 
tions for instantaneous and finite opening of diaphragms. 

Helium was selected as the test gas for the calculations so that perfect -gas rela- 
tions could be used without loss of accuracy. The initial conditions correspond to the 
data illustrated in figure 10(b), for which P4 is 32.7 MN/m^, T^ is 345 K, and pj is 
6.9 kN/m^ (run 585 in table II). To initiate the characteristic solution, the diaphragm 
was assumed to open in a stepwise manner; that is, the ratio of the cross-sectional area 
of the opened portion of the diaphragm to the tube cross-sectional area increased by 0.1 
in ten equally spaced time steps. (In other words, the diaphragm opening was modeled in 
a manner in which the open area increased linearly with time.) It was also assumed that 
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the partially opened diaphragm represented a Laval nozzle for which quasi-steady flow 
existed immediately upstream and downstream of the nozzle throat. The boundary con- 
ditions at the diaphragm location were determined from the upstream and downstream 
quasi -steady flow Mach numbers calculated by using the assumed time history of the 
area ratio. 

A diagram in the distance -time plane is shown in figure 23 for a diaphragm opening 
time of 200 psec. Each downstream -facing compression wave interacts with the inter- 
face and, upon reaching the shock wave, increases its strength. Thus, the shock- 
processed test gas displays a varying entropy because of varying shock strength with 
time. At the interaction of each compression wave with the shock wave, a reflected 
upstream -facing wave is produced. Although these relatively weak reflected waves were 
computed, their subsequent effect on the interface was ignored. In figure 23, the calcu- 
lated shock wave reaches a maximum velocity at a distance x of approximately 9 m. 
(The subsequent decay of the shock wave is not computed because of neglect of the 
reflection waves.) 

A finite -difference scheme to solve the one -dimensional Navier -Stokes equations 
for the finite diaphragm opening time problem was also used. This solution, generated 
by Weilmuenster (ref. 22), is a shock smearing method which uses a Lax-Wendroff differ- 
encing technique. The trajectory of the shock wave computed from this finite -difference 
scheme with an assumed opening time of 200 jusec is also shown in figure 23. As for the 
method -of -characteristics solution, the ratio of diaphragm opening cross-sectional area 
to tube cross-sectional area was assumed to increase linearly with time. The two theo- 
retical predictions of shock position as a function of time t are observed to be in good 
agreement. However, both predictions yield shock trajectories which trail experiment 
after the first 3 m or so. This is illustrated more clearly in figure 24. 

In figure 24, shock velocities predicted by the method of characteristics and the 
finite -difference scheme for helium as the test gas are compared with experiment. The 
characteristics solution and finite -difference solution yield values of I (location of the 
maximum shock velocity) somewhat larger than measurement. Neither of these solutions 
yield values of maximum shock velocity that exceed the predictions of references 10 
and 20, which are nearly identical for helium. Thus, all prediction methods considered 
significantly underpredict the experimental maximum shock velocities for helium. 
Assumption of diaphragm opening times greater than 200 psec in the characteristics 
solution and finite -difference solution will not alter the predicted maximum velocity. 
However, the location of maximum velocity will increase linearly with assumed opening 
time. 

Pressure -time histories of the expanding helium driver are shown at two axial loca- 
tions in figure 25. The measured pressure is higher at any given time than the pressure 
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predicted by conventional theory (instantaneous openii^ diaphragm) or by solutions which 
assume a 200-jusec diaphragm opening time. In fact, increasing the assumed opening 
time decreases the predicted pressure at any instance of time and thereby results in even 
poorer agreement with experiment. These pressure histories also support the conclu- 
sion that the finite opening time solutions do not adequately predict the present experi- 
mental results. 

Interface mixing. - White points out in reference 10 that a volume change occurs at 
the interface as a result of mixing of the driver and driven gases, if their specific heats 
are different. If the specific heat of the driven gas is greater than that of the driver gas, 
this volume change is positive and therefore accelerates the shock. Specifically, for a 
large temperature ratio across the interface. White derived the expression 


1 


V 21 



(5) 


where AV is the change in volume and c « and c r> are the specific heats at con- 
stant pressure of the driven gas and the driver gas, respectively. For air and CO 2 
driven by helium, excitation of internal degrees of freedom increases the effective spe- 
cific heat, and thus provides significant volume changes. However, for tests made with 
helium as the test gas, the helium flow behaves ideally, and although imperfect (inter - 
molecular force) effects are present in the helium driver, AV is essentially zero. The 
helium tests of figure 22(d) show that measured shock velocity exceeds predictions which 
take into account elevated driver temperature and formation from compression, even in 
the absence of appreciable volume change due to mixing. Although mixing may be a con- 
tributor to higher velocities than predicted for the present results, it is not a principal 
contributor. 

Nonequilibrium flow .- For low quiescent test-gas pressures, corresponding to low 
values of density in the post -shock region, the relaxation of dissociation and ionization 
may be slow compared with the particle residence time behind the shock. This results 
in a lower density ratio across the shock wave P 2 /P 1 (higher temperature in post-shock 
region) and thus a faster shock velocity for a given flow velocity behind the shock wave 
(conservation of mass considerations). The effect of nonequilibrium cannot be an impor- 
tant contributor to the difference observed in shock velocity between prediction and mea- 
surement because (1) shock velocities exceeding prediction are observed for values 
of Pj^ sufficiently high to suppress an extensive departure from equilibrium flow, and 
(2) the helium results, in which chemistry and ionization are not factors, also display 
higher measured shock velocities than predicted. 

Departures from one -dimensional flow .- In the previous theoretical analyses the 
shock -tube flow is assumed to be one -dimensional throughout. However, several factors. 


27 



particularly the diaphragm opening and the imsteady wall boundary layer, may cause 
significant departures from one -dimensional flow. Analysis of these effects on the flow 
has received only limited attention in the literature. 

When a scribed metal diaphragm first ruptures, the driver gas issues from the 
opened portion and expands as a free jet to fill the driven-section cross section. A num- 
ber of curved compression waves propagate downstream and coalesce to form a curved 
shock wave. The profile of this shock wave becomes less curved as it proceeds down 
the driven section. Such a shock formation process is discussed and illustrated in ref- 
erence 23. At the interface a mixing process is initiated by the jet -plume nature of the 
expanded driver gas. Near the diaphragm an oblique shock system is required to deflect 
the three-dimensional emerging jet to rectilinear flow. This shock system will exist 
until the diaphragm is fully open. In reference 24 the diaphragm opening process was 
assumed to be a circular aperture enlarging with time, and the resulting flow was treated 
two -dimensionally by a finite -difference technique. Both driver and driven gases were 
assumed to be ideal air, and diaphragm pressure ratios P 4 /P 1 of 10 , 100 , and 1000 
were considered. The results of this study were in good agreement with the shock 
strength predictions of reference 20 for pressure ratios of 10 and 100. At a pressure 
ratio of 1000, the two-dimensional shock strength prediction of reference 24 was slightly 
greater than the one -dimensional prediction of reference 20. Thus, the predictions of 
reference 24 demonstrate that the trend resulting from departure from one -dimensional 
flow is to increase the shock velocity for a given pressure ratio. 

The wall boundary layer represents another departure from one -dimensional flow, 
and examination of the interface velocities has suggested the possibility that mixing of 
the driver and driven gases by the boundary layer may be important. The velocity of the 
interface was not measured accurately in the present experiments; however, for one group 
of tests, an attempt was made to determine the interface velocity from the photomultiplier 
traces. The test gas was air and the driver pressure was about 54 MN/m^. Interface 
arrival was assumed to correspond to the rapid decay of the photomultiplier signal. The 
resultant determinations of shock-interface separation were in good agreement with the 
prediction method of reference 25. 

Average velocities of the interface and the corresponding measured static pres- 
sures are shown in figure 26 for the group of five tests. Also shown in the figure is the 
predicted variation of pressure with velocity for an isentropic unsteady expansion of the 
helium driver, in which the initial state was assumed to be p^ equal to 54 MN/m^ and 
T 4 equal to 330 K. Comparison of the theory with experiment shows that velocity is 
significantly underpredicted for a given pressure. It appears that some mechanism is 
present which imparts additional energy to the expanding driver gas. 
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The unsteady wall boundary layer is one possible source for this additional energy. 
In the shocked gas region ahead of the interface, the boundary layer entrains a significant 
portion of the mass of gas passing through the shock wave (ref. 25). This shock -heated 
gas tends to "surround" the driver gas as illustrated in the sketch below. This boundary 


Expansion region 



layer in the driver -gas region may in many cases be thick compared with the tube radius 
and is likely to be turbvilent. Extensive mixing between the cold (expanded) driver gas 
and the entrained driven gas, resulting in an appreciably warmer driver gas, may occvir. 
As this warmed driver gas expands, it will display a different pressure -Velocity relation- 
ship than that assumed for isentropic expansion of the pvire driven gas. 


The previous discussion is speculative and is not supported by calculation because 
of lack of an adequate flow model to represent the thick two-gas boimdary layer. No 
applicable studies have been found in the literature. Since the various one -dimensional 
analyses considered have not provided a totally adequate explanation for the present data, 
it appears likely that phenomena such as those discussed may play a contributing role. 


CONCLUDING REMARKS 


Because the incident shock velocity in the intermediate section of an expansion tube 
is a basic input to conventional expansion -tube theory for predicting test-section flow 
quantities, it is measured for each test. This measured shock velocity must be known 
accurately, since small uncertainties in velocity result in large vmcertainties in predicted 
test-section quantities. Existing methods for inferring shock velocity were studied by 
operating the Langley 6 -inch expansion tube as a shock tube with air, argon, carbon diox- 
ide, and helium as test gases. Unheated helium was employed as the driver gas, and 
most data were obtained at pressures of approximately 34 and 54 MN/m^. A range of 
pressure ratio across the primary diaphragm was obtained by varying the quiescent test- 
gas pressure from 0.0276 to 34.5 kN/m^. 

Shock velocities were inferred from microwave interferometer measurements, 
response of pressure transducers positioned along the driven section (time -of -arrival 
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transducers), and meastired tube -wall pressure. The microwave technique gives better 
definition of shock -wave characteristics (acceleration-deceleration trends) along the tube 
and is a viable method provided that microwave reflection occurs at the shock front. For 
the present range of conditions, the microwave proved to be rather limited in application. 
It is recommended that values obtained by the microwave method should always be veri- 
fied by one of the other methods. Shock velocities inferred from time -of -arrival trans- 
ducers and from measured wall pressures agree within the experimental accuracies. 

The time -of -arrival method avoids many of the shortcomings of the other two methods 
and is believed to be the most acceptable method of those examined. Measured shock 
velocities were observed to accelerate initially and then to decelerate with distance 
downstream from the diaphragm for all test gases. This trend was most pronounced for 
helium; argon, air, and carbon dioxide experienced successively less attenuation. The 
distance required for the shock velocity to reach a maximum value increased with 
increasing pressure ratio across the diaphragm for all test gases. At the higher pres- 
sure ratios, the shock velocity remained more nearly constant with distance after the 
initial acceleration for the present tube length of 29.0 m. For a given test gas and pres- 
sure ratio, minor differences in shock velocity were observed between single- and 
double -diaphragm modes of operation. The measured maximum shock velocity with 

O 

helium as the test gas and a pressure ratio of approximately 10"^ was 1.08 to 1.11 times 
the velocity predicted from conventional shock -tube theory for all diaphragm materials 
tested; however, the distance for the shock velocity to reach a maximum value varied 
with diaphragm material. 

Measured shock velocities exceeded real -gas predictions from conventional shock- 
tube theory and predictions from formation-from -compression theories for all test gases, 
with the exception of some argon results. Primary contributing factors to this discrep- 
ancy are believed to be (1) the effect of the helium driver temperature increasing upon 
pressurization of the driver section, (2) effects resulting from the finite opening time of 
the primary diaphragm, and (3) two-dimensional or nonplanar flow effects. Tests with 
helium as the test gas discounted interface mixing and flow chemistry as primary con- 
tributors to this difference between measured shock velocity and theory. 

Lar^ley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
November 20, 1975 
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TABLE I.- LOCATION OF TRANSDUCERS USED TO GENERATE 
TIME -OF -ARRIVAL DATA FOR BOTH SERIES OF TESTS 



^21.6 -m-long driven section. 
’^29. 0-m-long driven section. 
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TABLE n.- CALCULATED THERMOCHEMICAL-EQUILIBRIUM PROPERTIES BEHIND INCIDENT 
SHOCK FOR TESTS WITH DOUBLE DIAPHRAGM AND DRIVER PRESSURE OF 
APPROXIMATELY 33 MN/m^ FOR VARIOUS TEST GASES 


(a) Argon 


Run 

Tl, 

K 

Pi, 

kN/m2 

^s,max, 
km /sec 


At Ug 

,max 

Us, final’ 
km /sec 


At U 

s, final 

P2’ „ 
kN/m2 

T2, 

K 

3 

electrons/cm'^ 

P2’ ^ 

kN/m2 

T2, 

K 

electrons/cm^ 

' 571 

297.0 

34.5 

2.088 

1820 

4 186 

1.62 X loll 

1.807 

1360 

3201 

1.51 X 10® 

570 

297.6 

6.90 

2.743 

627 

7 030 

7.47 X 10 14 

2.271 

429 

4906 

2.02 X 10 12 

569 

297.0 

3.45 

3.094 

403 

8 720 

7.82 X lolS 

2.475 

256 

5777 

2.70 X 10 12 

566 

297.0 

.690 

3.523 

109 

9 986 

1.57 X lol® 

2.987 

75.1 

8135 

1.56 X 10 15 

567 

297.0 

.345 

3.688 

61.1 

10 160 

1.36 X 10 1® 

3.261 

45.9 

9068 

3.96 X 10l5 

593 

298.2 

.0690 

3.947 

14.4 

10 030 

5.70 X lol5 

3.758 

12.9 

9727 

4.10 X 10l5 

594 

297.0 

.0356 

4.069 

8.06 

, 9 939 

3.87 X lolS 

3.978 

7.66 

9817 

3.38 X 10l5 


(b) Air 

Run 

Tl, 

K 

Pi, 

kN/m2 

Us, max’ 
km/sec 


At Ug 

,max 

Hs,final’ 

km/sec 


At U 

s,final 

P2’ 

kN/m2 

T2’ 

K 

Ne, 

electrons/cm 2 

P2’ 

kN/m^ 

T2’ 

K 

Ne, 

electrons/cm2 

561 

297.0 

34.5 

2.057 

1470 

2060 

2.99 X lo"^ 

— 

1.905 

1250 

1833 

1.05 X 10® 

559 

297.2 

6.90 

2.606 

479 

2883 

6.73 X lolO 

2.332 

380 

2480 

2.17X 10® 

557 

297.8 

3.45 

2.911 

302 

3227 

4.63 X loll 

2.576 

233 

2819 

2.95 X lolO 

553 

298.7 

.690 

i 3.414 

84.5 

3543 

1.07 X 10l2 

3.100 

69.1 

3262 

2.52 X loll 

554 

297.6 

.345 

3.551 

46.1 

3573 

8.36 X loll 

3.283 

39.2 

I 3339 

2.70 X loll 

556 

298.2 

.0690 

4.100 

12.4 

4111 

1.94 X 10l2 

3.932 

11.4 

3807 

7.47 X loll 



Table n.- Concluded 


(c) CO2 



Tl, 

K 




At U 

s,max 



At U 

s, final 

Run 

kN/m2 

km/ sec 

P2. 

kN/m2 

T2. 

K 

Ne, 

electrons/cm^ 

^Sjiinai? 

km/sec 

P2» 

kN/m^ 

T2, 

K 

Ne, 

electrons/cm® 

580 

297.6 

6.90 

2.301 

598 

2285 

6.69 X 10^ 

2.149 

519 

2117 

2.96 X 10^ 

578 

297.6 

3.45 

2.530 

364 

2472 

7.91 X 10® 

2.347 

312 

2308 

7.66 X 10® 

573 

298.2 

.690 

3.155 

115 

2778 

1.72 X 10® 

2.819 

91.0 

2578 

1.79 X 10 '^ 

575 

297.6 

.345 

3.338 

64.6 

2814 

2.06 X 10® 

3.048 

53.6 

2658 

3.68 X lo"^ 

576 

297.0 

.0690 

3.804 

16.9 

2872 

2.16 X 10® 

3.642 

15.5 

2799 

9.75 X 10 '^ 

577 

297.0 

.0345 

4.084 

9.78 

2919 

2.25 X 10® 

3.962 

9.20 

2865 

1.37 X 10® 


(d) Helium 
Ug max 


Ug final 


Run 

Tl, 

K 

Pl> 

kN/m® 

^s,max, 

km/sec 

P2, ^ 

kN/m^ 

T2, 

K 

Ne, 

electrons/cm3 

11s, final, 
km/sec 

P2, ,, 

kN/m® 

T2, 

K 

Ne, „ 

electrons/cm'* 

586 

297.0 

34.5 

3.490 

502 

1355 

1.56 X 10-14 

2.865 

336 

994 

3.22 X 10-14 

585 

297.0 

6.90 

3.962 

130 

1673 

1.29 X 10-14 

3.307 

90.0 

1242 

5.00 X 10-1® 

582 

297.6 

3.45 

4.206 

73.2 

1854 

5.67 X 10-14 

3.597 

53.3 

1424 

3.67 X 10-1® 

583 

297.0 

.690 

4.511 

16.9 

2094 

1.48 X 10-10 

4.298 

15.3 

1924 

3.89 X 10-13 

584 

297.6 

.345 

4.511 

8.43 

2095 

1.06 X 10-10 

4.420 

8.09 

2021 

9.12 X 10-12 
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Figure 1.- Schematic diagram of expansion-tube flow sequence. 
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(b) Air. Single diaphragm; p^ = 0.690 kN/m^. 


to 

c 

(U 

CO 


I 




^ A .A / 
. J V 

V V ^ 


t ^,AA/ \ A r. A A A A «. . 




1 msec 


H 


(c) CO 2 . Double diaphragm; pj = 3.45 kN/m^. 
Figure 4.- Sample microwave traces for several test gases. 
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(a) Argon. 
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(b) Air. 
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(c) CO 2 . L-75-250 

Figure 5. - Sample oscilloscope traces of wall static pressure 
in driven section for several test gases, pj = 3.45 kN/m^; 

X = 19.75 m. 
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U^, km/ sec U^, km/ sec 

(a) Air. Ug = 2.8 km /sec. 


Figure 6.- Uncertainty in free-stream static pressure and pitot pressure due to corresponding uncertainty 
in shock velocity in intermediate section of expansion tube for several test gases, pj = 3.45 kN/m^; 

Tj = 300 K. Solid line was computed with program of reference 5 and shaded region denotes uncertainty 
corresponding to 5 percent uncertainty in Ug. 
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(a) P]^ = 34.5 kN/m^. 



X, tn 


(b) pj = 6.9 kN/m2. 
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(c) Pi = 3.45 kN/m2. 

Figure 7. - Incident shock velocity as a function of distance downstream 
of primary diaphragm for air. 
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(a) Pi = 34.5 kN/m2. 



(b) pj = 6.9 kN/m2. 



(c) Pi = 3.45 kN/m2. 

Figure 8. - Incident shock velocity as a function of distance downstream 
of primary diaphragm for argon. 
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(a) pj = 34.5 kN/m2. 
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(b) pj = 6.9 kN/m2. 



(c) Pi = 3.45 kN/m2. 

Figure 10. - Incident shock velocity as a function of distance downstream 
of primary diaphragm for helium. 
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(d) Pi = 0.690 kN/m^. 



X, m 


(e) Pi = 0.345 kN/m2. 
Figure 10.- Concluded. 
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(a) Air. 

Figure 11.- Time of shock arrival as a function of axial distance 
downstream of primary diaphragm. 
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(b) Argon. 

Figure 11.- Continued. 
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(c) COg. 

Figure 11.- Concluded. 










(a) Driver -gas pressure. 



(b) Driver -gas temperature. 

Figure 15.- Time history of helium driver pressure and temperature 
during pressurization of driver section. 




(a) Air. 

Figure 16.- Measured maximum incident shock Mach number in various test gases 
and prediction from conventional real -gas shock-tube theory. 
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Conventional theory 
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(b) Argon. 

Figure 16.- Continued. 
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Figure 16.- Continued. 
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(d) Helium. 

Figure 16.- Concluded. 
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(a) Air. 

Figure 18.- Measured incident shock velocities for double- and single -diaphragm 
modes of operation with various test gases. Measured Ug obtained from 
counter -timer readings; p^ « 34 MN/m^. 


67 


Double diaphragm 


Single diaphragm 



p-j = 0.069 kN/m 


■ 




p.| = 6.9 kN/m^ 


-J 1 1 i- 

15 20 25 30 

X, m 


(b) Argon. 

Figure 18. - Continued. 
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Figure 18.- Continued. 


69 

I 



km/sec 




X, m 
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Figure 18.- Concluded. 
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Figure 19.- Helium driver pressure as a function of time for single- and 
double -diaphragm modes of operation. 
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(a) Laminated plastic (polyethylene terephthalate). 



X, m 

(b) Stainless steel, AISI type 304. 


Figure 20. - Incident shock velocity as a function of distance downstream of 
diaphragm for various diaphragm materials. Helium as test gas; 

« 333 K; P 4 /P 1 “ 1030. Dashed lines denote conventional theory 
(ref. 5 ) with measured p^ and T 4 as inputs. 
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Figure 23.- Distance-time sketch of method-of -characteristics solution for helium 
driver and test gas. Diaphragm opening time of 200 /xsec. 
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Figure 24.- Experimental and theoretical variation of shock velocity with distance downstream 

of diaphragm for helium as driver and test gas. 
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Figure 26.- Measured and calculated interface conditions for a 
series of tests with air as test gas. 
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